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Abstract — We study the interference and outage statistics in a 
slotted Aloha ad hoc network, where the spatial distribution of 
nodes is non-stationary and isotropic. In such a network, outage 
probability and local throughput depend on both the particular 
location in the network and the shape of the spatial distribution. 
We derive in closed-form certain distributional properties of the 
interference that are important for analyzing wireless networks 
as a function of the location and the spatial shape. Our results 
focus on path loss exponents 2 and 4, the former case not 
being analyzable before due to the stationarity assumption of the 
spatial node distribution. We propose two metrics for measuring 
local throughput in non-stationary networks and discuss how our 
findings can be applied to both analysis and optimization. 

Index Terms — Ad hoc networks, interference, throughput 

I. Introduction 

The application of stochastic geometry to the modeling and 
analysis of wireless networks has attained a lot of attention 
in recent years. It has enabled a new framework called 
transmission capacity (TC) framework, which led to many new 
profound results in the topic of wireless networks (cf. JT], 0). 
The advantage of using a spatial model to describe the node 
positions rather than assuming a deterministic network topol- 
ogy is two-fold: First, such a probabilistic approach decouples 
the performance analysis from the actual topology, thereby in- 
creasing the generality of results. Second, it provides powerful 
means for network optimization, especially for highly dynamic 
networks, where interference is (unpredictably) fast-varying. 

With few exceptions, the node positions are mostly modeled 
as a stationary point process. Stationarity is a desirable prop- 
erty, allowing analytically tractable computations and, even 
more important, representing a key requirement for applying 
the definition of TC. Even though the stationarity assumption 
has not really narrowed the range of obtainable insights, it 
poses some shortcomings to the analysis of wireless networks: 
Infinite networks: Stationarity implies that the network is 
infinitely large as opposed to real deployments with a finite 
number of nodes. 

No border effects: Border effects are inherently neglected in 
infinite networks. However, border effects cause heterogeneity 
in the nodes' capabilities depending on their location, i.e., 
being disconnected, interference-/noise-limited, etc. 
Infinite interference for free-space path loss: For stationary 
node distributions in the plane and a path loss exponent a = 2, 
the interference is infinite almost sure (a.s.) (JT], resulting in 
a TC of zero. More specifically, stationary models lose their 
accuracy as the path loss exponent decreases due to the fact 



that infinitely many nodes contribute to the interference. 
Application of the TC: As already mentioned, the TC applies 
only to stationary networks. When the node distribution is non- 
stationary, this metric must be modified to take into account 
heterogeneous node deployments. 

In reality, wireless ad hoc networks always exhibit a het- 
erogeneous node distribution. The most obvious example is 
perhaps when the nodes are distributed in a bounded region. 
In such a network, the interference situation near the center 
will significantly differ from that at the border. Besides this 
simple example, more complex deployments are often found 
in practice, e.g., wireless sensor networks created by airdrop 
[3 1, spontaneous formation of hot spots [4|, etc. 

A. Contribution 

We extend the existing framework by relaxing the require- 
ments on the node distribution, thereby allowing for isotropy 
only. More specifically, we have the following results: 

• The interference and outage statistics for slotted Aloha 
with a = 2 and a = 4 are derived as a function of 
the receiver position and the spatial shape of the node 
distribution. We consider a path loss plus block fading 
model. As for the outage statistics, we focus on Rayleigh 
fading. We show how known results for the stationary 
case arise from our results as special cases. 

• Two global metrics, namely the differential TC and the 
average sum throughput, that take into account heteroge- 
neous node deployments are proposed. While the former 
metric is a refinement of the TC, the latter quantifies the 
first order overall network efficiency. 

B. Related work 

Stationary models with heterogeneous node deployment 
have already been investigated. Specifically, Poisson-Cluster 
and Matern hard-core models [6| have been studied, as 
they are well-suited for analyzing more sophisticated medium 
access control (MAC) schemes. Treated as general motion- 
invariant, these and similar models were further analyzed in 
ll2l. If7l. ll8l in a unifying way. In 191, a non-stationary and 
isotropic node distribution was assumed for analyzing multi- 
antenna receivers. While the analysis showed that the shape 
of the spatial distribution has a considerable impact on link 
performance, the scenario was limited only to the case of the 
receiver located in the origin. 



II. Network model 



III. Interference analysis 



We consider a wireless ad hoc network with nodes isotrop- 
ically distributed in R 2 . The MAC employed by the nodes is 
slotted Aloha. In a randomly chosen slot, some nodes wish to 
transmit a packet. We assume that the set of transmitters {x} 
follows an isotropic Poisson point process (PPP) <f> t := {x} 
on K 2 with intensity X(x), where x £ R 2 . Due to the isotropy 
of $ t , \(x) is rotation-invariant and depends only on the 
Euclidean norm ||x||, i.e., X(x) = A(||x||e^) = A(||x||), 
4> € [0, 2tt). When working with polar coordinates, we 
will use the notation A(r), where r := |jx||. With [6|, we 
can describe A(r) as the resulting intensity after distance- 
dependent thinning of a stationary PPP of intensity A, i.e., 



A(r) := AF(r), 



(1) 



where F(r) is called the shape function as it reflects the spatial 
shape of $ t . We will pose the following restrictions on F(r): 

(i) Positiveness: F(r) > for all r > 0. 

(ii) Normalization: max r {F(r)} = 1. 

The restrictions (i) and (ii) are necessary to ensure that A(r) 
is non-negative and bounded by A everywhere. 

We assume that each transmitter x has an intended receiver 
y randomly located at fixed distance d. From the random 
translation Theorem [6 1 it follows that the set of receivers {y} 
forms an isotropic PPP := {y} on R 2 with intensity X(x) 
as well. The fixed distance assumption is commonly accepted, 
see HI. However, we will relax this assumption in Section |V| 

We consider a path loss plus block fading channel with inde- 
pendent and identically distributed (i.i.d.) fading coefficients. 
The power path loss between two positions x, y £ R 2 is given 
by £(\\x — y\\) := (c + \\x — 2/|| Q )~ 1 with path loss exponent 
a. The parameter c > ensures boundedness of £. The power 
fading coefficient between a transmitter at x and a receiver at 
y is given by g xy , where E [g xy ] = 1 for all x, y £ R 2 . 

We further place a receiver at yo £ M 2 and an intended 
transmitter at an arbitrary position x$ £ R 2 with distance d 
to yo. The pair xo — > yo is called the reference pair as it will 
allow us to measure the (spatially-averaged) link performance 
for receivers at distance \\yo\\ from the origin. 

Assuming fixed power transmissions for all nodes, the 
instantaneous signal-to-interference-plus-noise ratio (SINR) at 
the reference receiver yo is given by: 



SINR(yo) := 



v +e(d)-i\( yo y 

where 77 is the average noise-to-signal ratio and 

xe$,\{z } 



(2) 



(3) 



is the interference power. We assume strong channel coding, 
i.e., the outage event is a steep function of the SINR. The 
outage probability (OP) at the reference pair xo — > yo is then 
given by the reduced Palm probability 

q(y ) :=P !a:o (SINR(yo) < P) , (4) 
where f3 is a modulation and coding specific SINR threshold. 



We now study the interference statistics at the reference 
receiver at yo- Before, we note the following two integral 
identities which are taken from iflOl : 

Identity 1. If a > \b\, a, b £ R, 



d<f> 



77 Pn ( v/^F?) 



(5) 



/„ (a + b cos ( a 2„ 5 2)i±± ' 

where P n {x) is the n th -Legendre polynomial. 

Identity 2. Let 01,02,03 £ R, R := a\ + a 2 t 2 + a^, A = 
4aia3 — a 2 , and 03 > 0. Using substitution t — > i 2 , we have 

' ' J A > (6) 

A = 0. 



V , == = <^ arcsinh 2 " 3 * + a2 , 
I log(2a 3 i 2 



a2), 

We are now in the position to derive the first moment of 
the interference at yo- 

Theorem 1. Let /(r) := dF(r)/dr, c > and a = 2. If 

lim F(r)r v < 00 for some v > 0, then 

1 — > 00 

E !a; "[l( 2/ o)] = AA 2 ( 2/ o,c)<oo, (7) 
where the interference-driving function A2(yo,c) is given by 



-^2(2/0, c ) : = ^(0) arcsinh 



f(r) arcsinh 



dr. 



(8) 



Proof: We write 
E te ° [l(y )] = A 



\£{\\x-y \\)F(x) dx 



what follows from Campbell's Theorem and Slivnyak's The- 
orem ifTTl . and from the i.i.d. property of g xyo . Applying 
Identity [T] and [2] yields the result. ■ 
The function A2(yo,c) in © has an interesting interpreta- 
tion: ^2(2/0, c) can be described as the interference field asso- 
ciated with the origin o, from which the remaining interference 
adds up differentially. 

Corollary 1. Summary of some special cases of Theorem Q] 
1) When we assume F(0) = 1 and f(r) < for all r £ R + , 
F(r) can be interpreted as a complementary cumulative 
distribution function (CDF) with respect to a random 
distance r to the origin, yielding 

-2 



M{yo,c) = arcsinh 



E 



arcsinh 



Vl 



2) Letting \\yo\\ — > 0, we further have 

A 2 (0, c) = log(l/2c) + E [log(2(r + c))] . 

3) Letting c — > 0, we have E [1(2/0)] = which is due to 
the resulting singularity of ^(||x — 2/0II) at x = yo, cf. J2). 



4) Sparse network (0 < lim r ^oo F(r)r y < oo, < v < 2): 
Remarkably, / °° rF(r)dr = oo but ¥/ x ° [\(yo)] < oo. 

5) Dense network (0 < lim^oo F{r)r v < oo, v —>■ 0): As 
expected 0, E [x ° [\(y )] = oo. 

1) has an interesting interpretation as well: The expectation 
can be seen as averaging the differential interference over 
r. Such an interpretation may be appropriate when analyzing 
networks with a priori unknown or fast-varying spatial config- 
urations, for which a CDF is then used to model their spatial 
shape. 4) implies \(yo) < oo a.s. although infinitely many 
nodes contribute to the interference on average. Note that 5) 
includes the homogeneous case with F(r) — 1 (/(?') = 0). 

We now extend the findings of Theorem Q] Before, we need 
the following Lemma. 

Lemma 1. Let a\,a% € R, 0,1 > 0. Then, 

2td(j)dt 



o a% + {t 2 + a 2 . — 2ta2 cos < 



2Ja{ 



= arctan 



2Re{n(t, 0,1,0,2)} 
1 - |re(i,ai,a 2 )| 2 ' 



where 



K(t,ai,a 2 ) 



- 1+J (t 2 +a 2 )) 2 



4t 2 a 2 



(9) 



(10) 



Proof: The basic idea is to decompose the integrand into 
partial fractions and to apply Identity Q] and [2] yielding ((9) after 
some algebraic manipulations. Note that according to [10], (H) 
and (O hold only for real-valued parameters. However, they 
were verified to hold also for complex-valued parameters. ■ 

Theorem 2. Let f(r) := dF(r)/dr, c > and a = 4. Then, 
E ,x °[\(y )] = \A 4 {y 0) c)<oc, (11) 

where A4 (y , c) is given by (TT2l below. 

Proof: The proof is analogous to the proof of Theorem Q] 
and uses the integral identity of Lemma Q] We further make 
use of (ii) in Section [XT] to show that linv^oo F(r) < 00. ■ 

Corollary 2. Summary of some special cases of Theorem [2] 

1) Case c — > 0: By taking the limit lim c _j.o ^4(2/0, c) in (fT2l . 
we observe that E- Xa [I (yo)] = 00, cf. 3) in Corollary |2] 

2) Homogeneous case: Let F(r) = 1. Then, f(r) = and 



lim arctan 



2Re{K(r,c,y )} 
1 - K(r,c,y )\ 2 



-f, a = 



00, 



yielding E- x ° [\(y )] = as expected, cf. 0. 

All results of this Corollary are consistent with the literature. 

The first moment of the interference is useful for bounding 
the interference distribution for the path loss only scenario. 



In case of Rayleigh fading channels, the Laplace transform 
of 1(2/0), i-e., A(j, )(s) := E [exp{-sl(y )}], is of significant 
importance, since it allows one to obtain the OP in closed- 
form. When treating the case a = 2, we will always assume 
that F(r) satisfies the additional condition of Theorem Q] 

Theorem 3. For g xy ~ Exp(l) for all x, y € M 2 (Rayleigh 
fading), the Laplace transform of l(j/o) at y E R 2 is 

£-Kyo)( s ) = exp{-As A a (y Q ,s + c)} , (13) 

for the cases a = 2 and a = 4, where ^(j/ojc) is given by 
© and A4(y , c) is given by dTzb . 

Proof: We write 



r l(w) ( fl ) ( ^E^ 



(6) 



n E &B „ [exP {-«fcwo^(l|x - Jtoll)}] 
_xe# t 



(l-£ g (^(||o;-j/o||))) A(x)(k , 



exp 



where (a) follows from algebraic manipulations and the i.i.d. 
property of the g xyo . (b) follows from the probability generat- 
ing functional and the Laplace functional of a PPP [6|. After 
noting that C g (s) = (1 + s) _1 for g ~ Exp(l), the integral is 
computed using Identity [T] and and Lemma [T] ■ 
Note that (a) in the proof holds for general point processes 
and some approximation techniques for computing the right- 
hand side already exist lfl2l . The (b) part is for PPPs only. 

Corollary 3. Setting F(r) = 1 for all r € M + and c = 0, we 
obtain the well-known result for the homogeneous case with 
a = 4 [2]: £ Kvo) (s) = exp{-\^-^}. 

IV. Outage and Local Throughput 

A. Outage probability 

We now study the OP for the reference pair xo — > j/o- In 
order to broadly discuss the impact of the spatial shape on 
the performance, we focus on the Rayleigh fading scenario. 
For other channel models, the interference moments derived 
in Section [Til] can be used to effectively bound the OP, e.g., 
using the Markov inequality [ 1] . We do not expect additional 
insights by considering also other channel models. 

Theorem 4. The OP for the Rayleigh fading scenario and a = 
2 respectively a — 4 is given by 



q(y ) = 1 - Ateo) (0(c + d a )) exp {-f3 v } 



(14) 



Proof: It is well-known that the OP for Aloha MAC and 
exponentially distributed power gains g xy can be written in 
terms of the Laplace transform of the interference [0, [0: 
We condition |2) on $ and evaluate the OP first with respect 
to g XoVo . We finally use (O with s — (3(c + d a ). ■ 
By means of (TPfl i in Theorem |4] we can now measure the 
OP for Rayleigh fading at an arbitrary location for an arbitrary 



Ai(y ,c) := ^-j= [F{r) arctan ^ 



2Rej>(r, c,y )} 



l«(r,c,2/o)|" 



r=0 



/(,) arctan C ' dr 

1 - |/c(r,c,2/o)r 



(12) 



spatial shape function F(r) satisfying the given restrictions. 
Fig. [T] shows 5OI2/0II) vs. \\yo\\ for a = 2 and a — 4, thereby 
confirming the analysis. It can further be observed how the 
network "moves" from the interference-limited to the power- 
limited regime with increasing \\yo\\. 

To highlight the accuracy of the model, we compare the OP 
from Theorem |4] to a straightforward way of approximating 
the OP which consists of assuming that the intensity X(x) is 
approximately constant around j/o- in this case the OP can then 
be described as in the homogeneous case [2|, except for the 
intensity in the exponential term being modulated by F(y ), 
i.e., q(y ) := 1 - exp{-F(y )Xir 2 d 2 pi § esc ^} « q(y ). 
We will now study the logarithmic ratio of exact to approxi- 
mate success probability, i.e., 7 := log 

Corollary 4. Let c = 0. The ratio 7 for a = 4 is given by 

7 = Xd 2 ^ (£F(yo) - d 2 ^A A (y a ,(3d 4 ) 



(15) 

Fig. [2] shows the ratio 7 together with the shape function 
F(r) for different receiver positions yo, F(r) was chosen such 
that the network exhibits a communication hotspot, with the 
density of active nodes slowly decaying between 70 and 500 
until it becomes approximately zero. One can see that the 
approximation is not satisfactory, especially in the transition 
region, where border effects come into play. 

B. Local throughput 

We now propose two local throughout metrics that are 
suitable for non-stationary wireless ad hoc networks. 

Definition 1 (Differential transmission capacity (DTC)). The 
DTC is defined as the maximal density of concurrent trans- 
missions in an infinitesimal region around the point x £ M 2 
subject to an OP constraint e, i.e., 



c(x,e) := X(x,e)(l - e). 



(16) 



The TC and its differential counterpart have similar mean- 
ing, except that the latter is position-dependent: For a given 
spatial shape F(r) and target OP e, c(x, e) yields the TC in a 
region dx. Hence, the DTC implicitly takes into account the 
spatial shape of the node distribution. For Rayleigh fading, 
c{x, e) is obtained by solving ( TT4i > for A. Like the TC, the DTC 
can be used for comparing different transmission protocols. 

Definition 2 (Average sum throughput (AST)). The AST is 
defined as the ratio of average number of successful transmis- 
sions to average number of simultaneous transmissions, i.e., 



E £ l{xs 



n := 



xG$, 



{x successful} 



E 



J2 l{x£R 2 } 
xG*, 



(17) 



The AST quantifies the first order overall efficiency of the 
network on the MAC layer. While the DTC highlights the 
spatial dynamics of the local throughput, the AST yields a 
single figure of merit. In essence, the AST counts the num- 
ber of successful transmissions, thereby integrating over the 
spatial dynamics. Note that the success function 1 { x successful} > 
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Fig. 1. q(y ) vs. \\y \\ for F(r) = exp{-(r/100) 3 }, d = 10, /3 = 0.5, 
c = 1, A = 0.001. Marks represent the simulation results. 
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Fig. 2. 7 and F(||l/o||) vs. \\y \\ for a = 4, p = 1, d = 10, c = 1. 

indicating that transmitter x has been successful, can be chosen 
arbitrarily to include additional outage-inducing effects, e.g., 
energy-limitations, dis-connectivity or secrecy outage. 

Theorem 5. Let lim r _ i . 00 F(r)r v < 00 for some v > 2. 
With the underlying network model and success function 
1{SINR(j/)>/3}> the AST can be computed as 

f °°r(l- q (r))F(r)dr 



fl = 



Jo rF(r)dr 



(18) 



Proof: Since the denominator directly follows from 
Campbell's Theorem, we focus on the numerator and write 



E 

(6) 
(c) 

(£) 

(e) 



^ , l{x successful} 
.x£$, 



(a) 



E'* [t{ x successful}] X(x)dx 

E$, [l{siNR(y)>«] P (y = y\x) dy X(x) dx 
E$, [l { siNR(y)>^}] P (y = y\x) X(x) dx) dy 
P| (SINR(y) >£)P(y = y\x) X{x)dx\ dy 

(l-?(l/))f/ F (y = v\z) Hx)dx) dy 
(l-q(y))X(y)dy. 



(a) is due to Campbell's Theorem ifTTI . (b) is obtained by 
noting that a transmitter x is successful if the intended receiver 
at y is not in outage. From Section [II] we know that y is placed 
by random translation of x according to some probability 



kernel P (y = y\x). (c) follows from Tonelli's Theorem Wi\ 
and (d) follows from E [l{xeA}] = F(X e A), (e) follows 
from (O and the fact that q(y) is independent of x, ■ 

V. Applications of the Model 

A. Shot-range inhibition 

Besides slotted Aloha, other MAC protocols such as 
CSMA/CA or local FDMA, are promising techniques for 
reducing excessive interference generated by nodes within 
shot-range. To study ad hoc networks with such inhibition 
mechanisms while ensuring analytical tractability, powerful 
methods based on non-homogeneous Poisson approximation 
have been used (6), fl4l . ITT5I . When such protocols are 
transmitter-initiated, e.g., transmitter sensing for CSMA or 
transmitter orthogonalization for FDMA, the resulting spa- 
tial distribution of interferers becomes inhomogeneous and 
approximately isotropic around the transmitter x, while the 
interference field at the intended receiver y will depend on 
the distance ||a; — y\\. Hence, our model can be applied also 
to such modeling problems and is not limited to Aloha MAC. 

B. Network optimization 

Consider the following situation: Let the set {y} of potential 
receivers be distributed as an isotropic PPP <5> r of intensity 
A r (r) = \F(r). Assume that <£> t and <& r are independent. That 
is, the set of all nodes follows a PPP, e.g., a sensor network 
created by airdrop, and connectivity at distance d is no longer 
guaranteed for every node. We further assume that the routing 
protocol employs a nearest neighbor strategy, i.e., transmitters 
aim at minimizing d. For points distributed as a PPP, the CDF 
i*d (d) of the distance d between a point and its nearest neigh- 
bor is well-known, see . We would like to know the optimal 
SINR threshold B such that the product log 2 (l + /3) Q(J3) with 
success function 1{ X successful} l{x connected} ls maximized. This 
corresponds to maximizing the expected sum rate, i.e., 

/T = argmax{log 2 (l+/?)ft(/3)} 

/■ \ f poo poo > 

« argmax|log 2 (l + /3)y rj (1 - q{r, 8,d))F d (dd)drl 

where we altered the notation f2 — ► fl(B) and q(r) —> q(r, /3, d) 
to point out the functional dependencies, (a) follows from 

E [9(l|y 4 ||)K = Xi] « q(\\x t \\), 

which essentially approximates the interference field at a 
receiver y^ by the interference field at the associated trans- 
mitter xi. This approximation is reasonable for high A r and/or 
moderate slopes of F(r). Fig.[3]shows log 2 (l + /3) f2(/3) vs. B. 
As can be seen, optimizing over B yields large improvements. 

VI. Concluding Remarks 

We extended prior work on the modeling and analysis of 
wireless networks by assuming an isotropic but not necessarily 
stationary spatial distribution of nodes. We derived, for slotted 
Aloha, the interference and outage statistics as a function 
of the receiver position and the shape of the spatial node 
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Fig. 3. log 2 (l + ft) f2(/3) vs. fj for a = 2. r)io denotes r\ at a distance 
d = 10. F(r) = exp{— r/250}. Marks represent simulation results. 

distribution. The case a = 2, which could not be studied yet 
due to the stationarity assumption, was intensively studied. 
For a = 4, we also obtain closed-form results, from which 
known results arise as special cases. We proposed two metrics 
for measuring local throughput in non-stationary and finite 
networks and discussed possible applications of our model. 
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